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Abstract 

In the present paper an influence of the anisotropic antisymmetric 
exchange interaction, the Dzialoshinskii-Moriya (DM) interaction, on en- 
tanglement of two qubits in various magnetic spin models, including the 
pure DM model and the most general XY Z model are studied. We find 
that the time evolution generated by DM interaction can implement the 
SWAP gate and discuss realistic quasi-one-dimensional magnets where it 
can be realized. It is shown that inclusion of the DM interaction to any 
Heisenberg model creates, when it does not exist, or strengthens, when 
it exists, the entanglement. We give physical explanation of these results 
by studying the ground state of the systems at T=0. Nonanalytic de- 
pendence of the concurrence on the DM interaction and its relation with 
quantum phase transition is indicated. Our results show that spin models 
with the DM coupling have some potential applications in quantum com- 
putations and the DM interaction could be an efhcient control parameter 
of entanglement. 



1 Introduction 

The entanglement property has been discussed at the early years of quantum 
mechanics as a specifical quantum mechanical nonlocal correlation [T]- [3] and 
recently it becomes a key point of the quantum information theory [4] . For en- 
tangled subsystems the whole state vector cannot be separated into a product of 
the subsystem states. This is why these subsystems are no longer independent, 
even if they are far separated spatially. A measurement on one subsystem not 



only gives information about the other subsystem, but also provides possibility 
of manipulating it. Therefore entanglement becomes the main tool in quantum 
computations and information processing, quantum cryptography, teleportation 
and etc. [5]. Due to the intrinsic pairwise character of the entanglement, entan- 
gled qubit pairs play crucial role in such computations. It is clear that single 
qubit gates are unable to generate entanglement in an N qubit system, and to 
prepare an entangled state one needs an inter qubit interaction, which is a two 
qubit gate. The simplest two qubit interaction is described by the Ising one be- 
tween spin 1/2 particles in the form of Jaf ctj. More general interaction between 
two qubits is given by the Heisenberg magnetic spin interaction models. These 
models have been extensively studied during several decades, experimentally in 
condensed matter systems [6] and theoretically as exactly solvable many body 
problems (Bethe, Baxter and others) [7], [8]. Now they become promising to 
realize quantum computation and information processing, by generating entan- 
gled qubits and constructing quantum gates [9], [10] in a more general context 
than the magnetic chains. 

Recently in this way interaction of two nuclear spins having the Heisenberg 
form were considered [11]. The nuclear spins from one side are well isolated 
from the environment and their decoherence time is sufRciently long. From an- 
other side nuclei with spin 1/2 are natural representatives of qubits in quantum 
information processing, which can realize quantum computational algorithms 
by using NMR [H], [H], [H]. 

Very recently entanglement of two qubits \l5i and its dependence on external 
magnetic fields, anisotropy and temperature have been considered in several 
Heisenberg models: the Ising model [l6],[T7], [H]; the XX and XY models 
[9], [19], [20], [21], [22], [23], [M]; the XXX model [25]; the XXZ model [26]; and 
the XYZ model [27] , [2H] , [59] ■ Particularly dependence of entanglement on the 
type of spin ordering, was shown, so that in the isotropic Heisenberg spin chain 
(the XXX model) spin states are unentangled in the ferromagnetic case J < 0, 
while for the antiferromagnetic case J > entanglement occurs for sufficiently 
small temperature T < = ki^s ■ Significant point in the study of such 
models is how to increase entanglement in situation when it already exists or to 
create entanglement in situation when it does not exist. Certainly this can be 
expected from a generalization of bilinear spin-spin interaction of the Heisenberg 
form. Around 50 years ago explaining weak ferromagnetism of antiferromagnetic 
crystals {a — Fe20^, MnCO^ and CrF^), has been controversial problem for a 
decade, Dzialoshinskii iJCF from phenomenological arguments, and Moriya [31j 
from microscopic grounds, have introduced anisotropic antisymmetric exchange 
interaction, the Dzialoshinskii-Moriya (DM) interaction, expressed by 

D-[SixS2]. (1) 

This interaction arising from extension of the Anderson superexchange interac- 
tion theory by including the spin orbit coupling effect [31, , is important not only 
for the weak ferromagnetism but also for the spin arrangement in antiferromag- 
nets of low symmetry. In contrast to the Heisenberg interaction which tends 
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to render neighbor spins parallel, the DM interaction has the effect of turning 
them perpendicular to one another. As we will see in the present paper it turns 
out that such spin arrangements are likely to increase entanglement. 

In most materials with weak ferromagnetism and the DM coupling, param- 
eter D is small compared to J. The values reported in the literature range 
from -J K, 0.02 to 0.07 (see [3^ and references therein). However in some com- 
pounds the DM interaction can attain a sizeable value in comparison with the 
usual symmetric superexchange J. Depending on compound its value varies be- 
tween -J « 0.05 to 0.2. Moreover, recently the DM interaction was found to be 
present in a number of quasi-onc-dimensional magnets [41] . Even it was found 
that the compound RbCoCl3.2H20 is described as a pure DM chain [5S]. The 
low-temperature magnetic behaviour of this compound gives strong evidence 
that the material consists of weakly interacting linear chains with predominant 
DM interaction. In addition, study of the DM interaction influence on dynam- 
ics of the one dimensional quantum antiferromagnet shows the big difference in 
the behaviour, depending on whether the coupling D is smaller or larger than 
the exchange interaction J [IT]. All these results imply that a study of spin 
models with DM interaction could have realistic applications. Then for applica- 
tions in quantum computations it poses the problem to find the entanglement 
dependence on this interaction. 

In the present paper we study the influence of the Dzialoshinskii-Moriya 
interaction on entanglement of two qubits in all particular magnetic spin mod- 
els, including the most general XYZ model. We find that in all cases, in- 
clusion of the DM interaction creates, when it docs not exist, or strength- 
ens, when it exists, entanglement. For example, we show that in the case of 
isotropic Heisenberg XXX model discussed above, inclusion of this term in- 
creases entanglement for antiferromagnetic case and for sufficiently strong cou- 
pling D > {kTsinh~^e^'^^^''^ — J^)^/^ it creates entanglement even in ferromag- 
netic case. We give detailed physical explanations of these results by studying 
ground state of the system at T=0. In this state we find nonanalytic dependence 
of concurrence on the DM interaction and establish its relation with the quan- 
tum phase transition. In addition, we show that time evolution generated by 
DM interaction can be implemented as the SWAP gate. These results indicate 
that spin models with DM coupling have some potential applications in quan- 
tum computations, and DM interaction could be an efficient control parameter 
of entanglement. 

The paper is organized as follows. In Section 2 we formulate the general 
XYZ model with DM coupling and find the density matrix and eigenvalues 
for the concurrence. Then we consider the time evolution and its relation with 
the SWAP gate. Since the concurrence calculation depends on several param- 
eters, in the following sections we consider all possible particular cases from 
the unified point of view. We think that such presentation is pedagogical and 
could be affordable by experimentalists. In Section 3, the main properties and 
entanglement of pure DM model and the relation of this model with SWAP 
gate are considered. The Ising model with DM interaction is studied in Section 
4. In particular, realization of the model for description of two nuclear spins 
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with DM coupling and implications for the quantum phase transitions in the 
presence of magnetic field are given. In Section 5 we consider the XY model 
and its particular reductions to the XX case, and to the Ising model. We show 
that inclusion of the transverse magnetic field leads to the different behaviour of 
concurrence C12 for the undercritical and the overcritical couplings. For T = 
the nonanalytic behaviour for C12 (D) is found. The XXX Heisenberg model is 
subject of Section 6. Section 7 is devoted to the XXZ model, where the influ- 
ence of DM coupling and magnetic field on the concurrence and the quantum 
phase transitions are studied. In Section 8 we study XYZ model in both anti- 
ferromagnetic and ferromagnetic cases, with inclusion of the DM coupling. The 
nonanalytic behaviour at T = is found. In Conclusions several implications 
for future studies are discussed. 

2 XYZ Heisenberg Model 

We start our consideration with the most general XYZ model, by inclusion of 
homogeneous B and nonhomogeneous b magnetic fields, and choosing the DM 
interaction ([T]) in the form -y = -j • Then for two qubits we have Hamiltonian 

^[J^ a^a^ + Jy a\al + J, ajd^ + B+ + B_ + Diafal - a\a-^)] (2) 



H 



where B+ = B + b, £?_ = B — b and erf , erf , cr|, i ~ 1,2 denote Pauli matrices 
related with the first and the second qubits. 



2.1 Eigenvalues and Eigenvectors 

To study the thermal entanglement in this system firstly we need to obtain 
all eigenvalues and eigenstates of the Hamiltonian H\'i>i) = Ei\'^i), i ^ 
1, 2, 3, 4. Simple calculations show that the energy levels are: 



7-^ '^Z 7-^ Jz 



(3) 



where = ^ B"^ ^ Jl,v = ^b"^ + J'l + D^, J± 
wave functions are 

J_ 

1 



-, and the corresponding 



1*1,2) 








1*3,4) 



VMv T b) 





(bTv) 
J+ - iD 




(4) 



For B = Q, b = 0, D^O these wave functions reduce to the maximally 
entangled Bell states 



1*2,1) — iBo,3) = ^(|00)±|11)) 
1*4,3) ^ |Bi,2) = -i=(|01)±|10) 



(5) 
(6) 
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2.2 Time Evolution of States and SWAP Gate 

Here we like to show the direct relationship between our spin model and quan- 
tum gates. For this we consider the evolution operator 



U{t) = exp[--Ht] 



(7) 



determined by two qubit Hamiltonian ^ of XYZ model with DM coupling, 
i? = 0, 6 = 0. Then evolution of the standard basis is given by 



100) 

111) 
101) 
|10) 



e 2R 



iJzt 

e 2R 



cos^lOO) -zsin%|ll) 



cos^|ll)-isin^|00) 



tv , , J+ — iD tv , 

cos— 101) -i-t sin— 110) 

n v h 

tu , , Ja- + iD tv , 

cos— 110) -i-^— sin— 101) 

n V h 



(8) 
(9) 
(10) 

(11) 



where v = ^ J'^ + D'^. In particular cases, discussed in the next section, this 
evolution can implement the SWAP gate at time t — fnr/2v. 



2.3 Density Matrix and Concurrence 

State of the system at thermal equilibrium is determined by the density matrix 

„-H/kT p-H/kT 



P{T) 



(12) 



where Z = Tr[e ^/'^^j is the partition function, k is the Boltzmann constant 
and T is the temperature. Then by exponentiation of Hamiltonian ([2]) we find 



-H/kT 



An 








Al4 





A22 


A2i 








As2 


A33 





A41 








Am 



(13) 



where 



All 

A44 

Ai4 

A41 



= e ^fci- 



= e 2fcT 



2kT 

^ J- 



A* 


sinh 


kT " 


fi kT 




B . , fJ. 
H smh — 


kT 


/i kT 


sinh 






kT 


sinh 






kT 



(14) 
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and 



A22 

A32 
Z = Tr\e 



e2kT 



cosh 



V 


b 




V 






sinh 




kf ^ 


V 




V 


b 


sinh 


V 










V 




kf 



;2feT gmh ■ 



iD 



sinh 



kT 

V 

kf 



-H/kT 



1 = 2 [e cosh + e 2fcT cosh -^1 
'I kT kTi 



(15) 



(16) 



As p(T) represents a thermal state, the entanglement in this state is called the thermal 
entanglement. The degree of entanglement could be characterized by the concurrence 
C12, which is defined as [15], [32] 



C12 — max{\i — A2 — A3 — A4, 0}, 



(17) 



where Ai > A2 > A3 > A4 > are the ordered square roots of eigenvalues of the 
operator 

P12 =p(a^®a^)p*((T«®T''). (18) 

The concurrence is bounded function < C12 < 1, so that when C12 ~ 0, the states 
are unentangled, while for C12 = 1, the states are maximally entangled. 
For the general Hamiltonian ([2]) we find : 



Al,2 



A3, 4 



—— W iH smh"' T—T — smh -f- 

Z M fi^ kT^ fi kT 



e2kT 



1 + — i^— smh^ — 

kT 



^ u ""^fcT 



(19) 



Then, to calculate the concurrence we need to order these eigenvalues. Since they 
depend on several parameters, before studying the most general case, it is useful to 
treat all particular cases separately to clarify the influence of the DM coupling on the 
entanglement. Starting from pure DM model we study various Heisenberg models, 
including the general XYZ case. 

Before this, we like just to stress here the general observation on the concurrence 
H17[). If the biggest eigenvalue say Ai is degenerate, then its positive contribution 
would be compensated by the another degenerate one, so that C12 = and states 
are always unentangled. We will encounter this situation in several cases and it has a 
simple physical explanation. The degenerate biggest eigenvalues of the density matrix 
correspond to the minimal values of the energy, so that the ground state of the system 
becomes degenerate and no entanglement occurs. 



3 Pure DM Model 

3.0.1 Main Characteristics of DM Model 

As we discussed in introduction some realistic quasi-one dimensional compounds with 
predominance of DM interaction can be described as a pure DM chain :39 . Here we 
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consider the main characteristic properties of the DM couphng between two qubits 
and its influence on the entanglement. If in Hamiltonian Q we put Jx = Jy = Jz = 
and B — b — then the model is determined completely by the DM term ((T}. In 
this case the first two eigenstates become degenerate Ei = E2 — and E^^i = ±D. 
For definiteness we choose D > 0, then for T — the ground state of the system 
with energy E4 = —D is an entangled state |10) — i|01). When temperature increases 
this state becomes mixed with the higher states and entanglement decreases. But for 
sufficiently large value of D the ground state can be alienated so that entanglement 
increases. This shows that for a given D there exists kTc — -D/ln(l + \/2) so that 
for the under critical case T <Tc the states become entangled and the concurrence is 

C12 = ''p — ^ (See Fig.l). For T = the concurrence C12 ~ 1 and the ground state 
is maximally entangled. 



C12 




kT 



0.5 1 1.5 2 



Figure 1: Concurrence versus temperature for D ~ \ and Tc — 1.136 



3.0.2 DM Model and SWAP Gate 

The time evolution in pure DM model from one side is related with the SWAP gate, 
from another side can create maximally entangled states. In this case according to 
((8)l-(|ll[l for time evolution we have 

t/(g)|00) = 100), t/(^)lll) = lll) (20) 

U{^m) = -|10), t/(^)|10)=-|01). (21) 

Therefore we can see that the operator t/(f^) acts as the SWAP gate. Moreover at 
time t = -nh/AD the states |01) and |10) becomes maximally entangled Bell states. 

U{^)\Ql) = -L(|oi)-|io)) (22) 
t/(^)|10) = -L(|io) + |01)) (23) 

4 Ising Model 

For Ja: = Jj, = 0, Jz 7^ and B = h = Q,D = the Hamiltonian ((2} describes the Ising 
model . It was observed before that for pure Ising model in both the antiferromagnetic 
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,Jz > 0, and the ferromagnetic cases, Jz < 0, the concurrence is zero and the states 
are always unentangled [TB], [T7], |18] . The physical insight of such behavior is easy 
to understand. When J_ = J+ = the density matrix p H12|) is diagonal in the 
standard basis which implies the absence of quantum correlations. Despite of having 
four maximally entangled states as the eigenvectors, the states |\E'i,2) and |'I'3,4) are 
degenerated, so that the Ising thermal state has no entanglement. The situation 
does not change if one includes homogeneous B or nonhomogeneous b magnetic fields, 
because the density matrix p is still diagonal and no entanglement occurs . 



4.1 Ising Model with DM Coupling (5 = 0, 6 = 0, D ^ 0) 

In contrast to magnetic fields, which does not create entanglement, inclusion of the 
DM coupling contributes to the nondiagonal elements of p and creates entanglement. 



4.1.1 Antiferromagnetic Case (J^ > 0) 

In this case addition of the DM coupling to the Ising model splits the degenerate ground 
state with E3 — E4 = so that it becomes a singlet with E3 = — — D, for 

D > or £4 = — + -D, for _D < 0. At T = this leads to the maximally entangled 
state with C12 ~ 1. When temperature increases the maximally entangled ground 
state becomes mixed with the higher eigenstates and the entanglement decreases. 
However, for a given temperature by increasing the coupling D > Dc, where Dc = 
fcTsinh"^ e^'^^'^*'^, we can decrease this mixture and increase entanglement, so that 
the concurrence is 

sinhM„e-^.A^ 
C12 = ^1 — . (24) 

cosh ^ + 6"-^^/'='^ 

4.1.2 Ferromagnetic Case (J^ < 0) 

In this case the ground state for small D at T = is also a doublet and no entanglement 
occurs. However, with growing D the eigenstate E3 — — D is lowering so that at 
critical value Dc ~ \Jz\ the ground state becomes triplet. When D > Dc the ground 
state i?3 is maximally entangled singlet. With growing temperature, a mixture of 
this state with the higher states decreases entanglement. For given temperature T, 
there exist the critical value Dc = \ Jz\ + ^ ln(l + e"^'-^^ 1/'''^) so that for D > Dc the 
concurrence is 

sinhM_el-/.IAT 
C12 = — — . (25) 

cosh ^ + ei-'^i''''^ 

Comparison of (|24|l and (|25p shows that in the antiferromagnetic case the states 
can be entangled more easily than in the ferromagnetic one. 



4.2 Ising Model for Two Nuclear Spins with DM Coupling 

As an application of the above calculations here we discuss entanglement of two nuclear 
spins. Recently two nuclear spins were considered in a model with weak Heisenberg 
type interaction in a constant longitudinal magnetic field along z direction [TT] 

H = Hz+Hxy (26) 
Hz = — -(tJiCTi + a;20-2 + Jo-iCTa) (27) 
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Hxy = -=^{JoW2+Jo\al) (28) 

where the isotropic form for the spin coupling J is assumed, and loi,2 = (B ^ b) are 
the Larmor frequencies of two nuclear spins, h — 1. In the experiments, two different 
nuclear spins are selected, u>i 7^ UJ2 (we assume wi > UJ2), and the longitudinal con- 
stant magnetic field is in the order of IT Hz, so that aii,a;2 are much larger than J 
and 77 = ^^^^^^^ ^ 1. Hxy is non-diagonal in representation and due to quantum 
fluctuations of order rj^ , can be ignored. Thus, the Ising part Hz of the Hamiltonian is 
a well precise approximation [11]. However as we have seen above, for the Ising model 
with external magnetic fields no entanglement occurs, this is why two nuclear spins in 
this model are unentangled for any uo\ and U2- From another side, as follows from our 
consideration in Sec. 4.1 the addition of an interaction between qubits in the form of 
the DM coupling could make them entangled. Now by adding the DM interaction to 
two nuclear spin Hamiltonian ((27} we get the Ising model with homogeneous magnetic 
field B, nonhomogeneous magnetic field b and the DM interaction D. In the antifer- 
romagnetic and the ferromagnetic cases, when Jz ~ ±i Jz| respectively, for sufficiently 

strong D > Dc, where , sinh -^ ^'L^'' — e^^rf^, the states become entangled 

and the concurrence is 

^ sinh - e^^-ir 

C12 = jj^ (29) 

cosh -I- cosh T^e^ 



where B = {uji + tJ2)/2, 6 = (wi - lj2)/2 and v = ^i)'^ ^ ]j2 is worth 

to note that the homogeneous magnetic field B does not change critical value for 
the entanglement, but could change level of the entanglement. Moreover, increasing 
magnetic field, decreases value of the entanglement. It turns out that for the system 
at T = 0, the concurrence becomes nonanalytic when D — Dc 



C12 



u = Bt\M; (30) 

0, iy<BT\J.\, 



D_ 

2v ' 



that implies quantum phase transitions at the critical value Dc = (B \ Jz 



\'2 



4.2.1 Ising Model with DM Coupling and SWAP Gate 

li Jx = Jy = but Jz and D are nonvanishing and related by Jz = 8nD, (n = 
±1,±2...), then again like in Sec. 3. 02 the evolution operator U(jih/2D) acts as the 
SWAP gate. Our consideration shows that the Ising model, which was derived in 
several physical situations for interaction of qubits, with addition of the DM coupling, 
from one side leads to entanglement of states, from another side it can model the 
SWAP gate (f^ . (PT|) . This result shows that the Ising model with DM coupling have 
some potential applications in quantum computations. 



5 XY Heisenberg Model 

In the pure XY Heisenberg Model Jz = 0, / Jy and B = 0, 6 = 0, D = in (HI, for 
the antiferromagnetic case > 0, Jy > the ordered eigenvalues are A3 > Ai > A2 > 
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sinh cosh 

A4 and for sinh ^ > cosh the entanglement occurs with C12 ~ ■j^ 

cosh -^^+cosh 

In the ferromagnetic case Jx < 0, Jy < the entanglement occurs when sinh > 



cosh^ with the concurrence [23], [21], [24] 



sinh — cosh ■, „ , ^ 

C12 = f^. (31) 

cosh + cosh -j^ 



For the particular case of pure XX model, when Jx — Jy = J, in both antiferro- 
magnetic and ferromagnetic cases the states become entangled at sufficiently small 
temperature 

T<n^ (32) 
ksmh '^1 

As was shown in [9], |10].|19j. [27) inclusion of the magnetic field does not change this 
critical temperature. 

5.1 XY Heisenberg Model with DM Couphng (S = 0, 6 = 

0,D^0) 

By addition of the DM coupling eigenvalues become Ai, 2 ~ - — , -^3,4 ~ z 



where Z = 2 



I j_ I ,/ ,72+132 

cosh ^TT^i + cosh ^ . „ — 



kT ' kT 

In the antiferromagnetic case for any temperature T we can adjust sufficiently 



J? +-D2 



strong DM coupling D so that for sinh ^ — > cosh the entanglement occurs 
with concurrence 

sinh ^ t„ — cosh ^ , , 

C7i2 = ^. (33 

cosh -^^ — jj7j5 h cosh -pjT 

The ferromagnetic case gives the same result as the antiferromagnetic one. The com- 
parison with pure XY model shows that the level of entanglement increases with 
growing coupling D. 

In particular case Jx = Jy = J, the ordered eigenvalues are A4 = ^-^ — > A3 = 

^^^^r— > '^^1.2 = where = y^J^ + and Z = 2{l + cosh^). Then the 

sinh -jpf — 1 



entanglement occurs when sinh > 1 and the concurrence is C12 = , ,, 

"-^ cosh 1^ + 1 

Comparison with the pure XX model (|32p shows that the critical temperature 



k smh '^1 

in this case increases with growing D. For D — 1^3) in Q is the ground state 
with eigenvalue -E3 — —\J+\, which is maximally entangled Bell state, so that the 
concurrence C12 = 1. As T increases the concurrence decreases due to the mixing of 
other states with this maximally entangled on^ . 



^In ref. |26) entanglement in XX model with DM coupling was derived but not in the 
general XXZ case as it is claimed in the paper. 
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5.2 Ising Model in Transverse Magnetic Field 

As a particular case of the general XY model now we consider the transverse Ising 
model, when Jy = 0, with external magnetic field B in z— direction [21], and with 
addition of DM interaction: 



(35) 



The corresponding eigenvalues and the partition function Z can be written as follows 

+ J2 



A3. 



J2 



sinh^ 

kT 



J 

VB^ + J2 



sinh - 



kT 



1 T ^''^^°^ 



kT 



cosh tt:: h cosh 



(36) 



(37) 



(38) 



kT kT 

To find the maximal eigenvalue we compare the difference of A4 and A2 as a function 
of B,D and T, A4 - A2 f(B,D,T): 



J2 



■ sinh' 



+ J2 



J 



When /(B, D, T) = we find the critical D = Dc{B,T) as 
Bc(i3,T) = 



. , Vi?2 + J2 
smh ^j;^^ (39) 



(40) 



-J2 +T2 In 



1 + 



J2 



£2 + J2 



■ sinh' 



VB2 ^ J2 
kT 



J 



VB^ + J2 



: sinh 



VB2 ^ J2 
kT 



In Fig. 2 we plot Dc as a function of T for different values of magnetic field B = 
0.05, 0.5, 0.7, 1 (J = 1, fc = 1). The 3D plot of as a function of B and T for the 
same values of parameters is given in Fig. 3. 




Figure 2: Dc versus T for B = 0.05, 0.5, 0.7, 1 

For critical D = Dc, the eigenvalues are degenerate A2 = A4 and as a result the 
concurrence Ci2{B, Dc,T) — 0. However the value of concurrence is different for the 
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Figure 3: 3D plot Dc versus B and T 



under critical and the over critical cases. In under critical case when D < Dc the 
maximal eigenvalue is A2 and for the concurrence we have 

sinh - cosh 

C12 = max{ == ,0}, (41) 

cosh ^ h cosh ^ — 

while in the over critical case, when D > Dc, the maximum eigenvalue is A4 and the 
concurrence is 



smh W 1 + smh ^ 

= maa;{ ^=1= -j== , 0} . (42) 

cosh ^j, ' h cosh ^j, ' — 

In pure Ising model when _B = and D = as we can see from (|39p we have /(O, 0, T) = 
and no entanglement occurs. But as reported in [21) an addition of the transverse 
magnetic field to the Ising model could create entanglement. Now we can generalize 
these results by analyzing in addition the influence of DM interaction on entanglement 
in the Ising model with the magnetic field. When B = Q the addition of solely DM 
term creates entanglement at sufficiently strong D, and this value of D becomes bigger 
for higher temperatures. If we have both terms B ^ and D 7^ 0, then with increasing 
D the behavior of entanglement becomes nontrivial. In Figs. 4. a, 4.b, 4.c we show 
behavior of entanglement as a function of D for different temperatures. When T — 
entanglement is nonanalytic function of D, given by the step function 

D <Dc; 



Ci2{D) = { 0, D = Dc; (43) 

1, D> Dc . 

where Dc = B (see Fig. 4-a). This nonanalytic behavior signals on the quantum 
phase transition (35] appearing at D = Dc = 1. In Fig. 4-b at temperature T = 0.5 
the entanglement as a function of D decreases down to zero and at Dc ~ 0.75 reaches 
its nondifferentiable minima. After this it increases monotonically with growing D. 
For higher temperature T = 1 in Fig. 4-c, the entanglement is zero until D becomes 
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sufficiently strong ai D — Dc, where entanglement appears and monotonically grows 
with growing D. 




Figure 4: Concurrence of Ising model in transverse magnetic field versus D, 
when S = 1 and T = 0.01, 0.5, 1 



6 XXX Heisenberg Model 

In pure XXX model Jx — Jy = Jz = J and B — b = D — in entanglement 
behavior for the ferromagnetic and the antiferromagnetic cases is different. In the 
spectrum of the model we have three degenerate eigenstates with eigenvalue J/2 and 
one eigenstate with eigenvalue — 3J/2. It was observed before [55] that for the ferro- 
magnetic case ( J < 0) the concurrence is zero and the states are always unentangled. 
It happens because when J < 0, the ground state of the system is an equal mixture of 
the triplet states with energy, Ei = E2 = E4, — . The density matrix p is diagonal 
and inclusion of magnetic field does not change the result. Increasing temperature T 
just increases the singlet mixture with the triplet, which can only decrease entangle- 
ment [25] . [34] . The situation is different for the antiferromagnetic case when J > 0. 
In this case the ground state is the maximally entangled singlet state with E3 = 
so that the concurrence C12 = 1 at T = 0. It decreases with T due to mixing of the 
triplet higher states with the singlet ground state. For a given coupling constant J 
entanglement occurs at temperature T < [26| . 

6.1 XXX Heisenberg IModel with BM Coupling {B = 0,b = 

Now by adding DM coupling for the antiferromagnetic and the ferromagnetic cases, 
for J = ±|J| respectively, for a given temperature T the entanglement occurs when 
D > Dc = \/ (fcTsinh"""" eTl''l/'='^)2 — with the concurrence 

sinh^^-e^l^l/'^^ , , 

C7i2 = . (44) 

As we can see inclusion of the DM coupling, in the XXX model, increases entangle- 
ment in the antiferromagnetic case and creates entanglement even in the ferromagnetic 
case. This can be explained if we consider the eigenvalues of our Hamiltonian varying 
with D. 

For the antiferromagnetic case the ground state of the system remains singlet with 
energy £3 = — ^ — \/ J'^ + D"^ , while from degenerate excited triplet state one of the 
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energy levels _B4 = — + y/ is splitting up. With increasing coupling D the 

gap between ground state and the first excited doublet state is increasing, this is why 
the system becomes more entangled. 

In the ferromagnetic case, from unentangled triplet ground state one of the states 
splits with the energy £3 = ^ — V + D"^ ■ Then at temperature zero this state 
becomes maximally entangled ground state. This way the DM interaction creates 
entanglement in the ferromagnetic case. With increasing D the gap between singlet 
ground state and the first doublet state increases, this is why entanglement in the 
ferromagnetic case increases. 



7 XXZ Heisenberg Model 

When Jx = Jy = J 1^ Jz the Hamiltonian ((2]) becomes 

H = 7:[Jicri(^2 + o-io-f + A 0-^0-2) + B+cTi + B- al + Dialal - o-j'tJa)]. 



(45) 



where A = Jzj J ■ 

• In a pure XXZ ferromagnetic model when < and — | Jz| < J < j J^j or j A| > 
1, we have the degenerate maximal eigenvalues Ai = A2 and no entanglement 
occurs. This happens since the ground state of the system is doublet with 
eigenvalues Ei = E2 — 

• In particular case |A| = 1 or \J\ = \Jz\ we have reduction to the XXX model, 
where the energy level E3 merges to the ground state, and the last one becomes 
triplet state, as we discussed above in Sec. 6. 

• For J > and A > — 1 the maximal eigenvalue is A3 and the states are entangled 



when sinh -pf > e -^^Z*^ with the concurrence 
^ sinh 



cosh^ +e-J-/'=T' 



(46) 



• For J < and A < 1 the maximal eigenvalue is A4 and the states are entangled 

(47) 



for sinh \^ > e ^^/^^ with the concurrence 



sinh 


J\ 

kT 


_ g~J./kT 


cosh 


J\ 
kT 


^-J^/kT 



7.1 XXZ Heisenberg Model with DM Coupling (B = 0,b 

With addition of the DM coupling we have the eigenvalues 



Al,2 



l + e././feT(,osij 



A 



,^^J^+D^/kT 



3,4 



g-J^/kT gosh 



(48) 



Then for Jz < and | J^] > |J|, there exists critical value Dc = \/ — J"^ so that for 

he state 



D > Dc and sinh ^ ' > e ■'^/'^^ the states are entangled with the concurrence 



C12 ~ 



sinh ■ 



J^I/fcT 



cosh^^±^+elJ^I/'=r 



(49) 
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This happens because for Jz < 0, |Jz| > |J| and D — 0, the ground state is doublet 
with El = E2 = "■^^i ^iid by increasing D so that D > Dc, the higher energy level 

lowers to the singlet ground state which is maximally entangled. Comparison of 
H49p with (|47p shows that with growing D entanglement increases. 

It is worth to note that the concurrence 1)49^ for both signs of J is the same. 
Moreover, as easy to see in (I49|) parameters J and D appear symmetrically. It means 
that the concurrence could be increased by growing J with fixed D either by growing D 
with fixed J. This refiects the known result [35] on equivalence of the Heisenberg XXZ 
model with DM coupling to pure XXZ model with modified anisotropy parameter and 
a certain type of boundary conditions. In fact comparing entanglement in our formulas 
for pure antiferromagnetic case (|47[) with the one including the DM interaction (I49p . 



we can see that the concurrences are connected by the replacement J ^ Jy 1 + -j^-, 
which corresponds to the substitution for the anisotropy parameter in the pure XXZ 
model as A —=^==. 



7.2 XXZ Heisenberg Model with DM Coupling and Mag- 
netic Field 

If we take into account the DM interaction D and magnetic field B simultaneously, 
the above results for critical value of the DM coupling are still valid, but the level of 
entanglement decreases according to 



sinh^^5P^-e-'^^/'=^ , , 

C12 = . (50) 

cosh ^^^i^ + e-J'/^T cosh ^ 



For T = Q and > we have nonanalytic behavior 

f 1, y/TPTJ^ > B 



Ci2^{ |, VWTJ^ = B~J,; (51) 



0, VW+J^ <B-.h. 

which signals appearance of quantum phase transitions. The concurrence versus tem- 
perature for different values of coupling D is shown in Fig. 5, where J = 1 , Jz = 0.5 
and magnetic field B = 2. As we can see in general the entanglement decreases with 
growing temperature. However we like to emphasize that for D < Dc in Fig. 5a, when 
D = 0.1, the entanglement is increasing with growing temperature. This phenomena 
can be explained by the fact that for such values of the parameters at T = the ground 
state is the separable state with energy Ei — ^ — B = —1.75, and the concurrence is 
zero (see the last case in eqn. ((51}). When temperature increases the entangled state 
with energy E3 = =p V + = —1.255 becomes involved into the mixture and 
entanglement is increasing. 

When D — Dc the entanglement decreases smoothly from C12 = 0.5 (Fig. 5b, 
Dc = 1.118). By increasing D {D = 1.19), first it gives sharp decrease from C12 ~ 1 
(Fig. 5c) and then it vanishes slowly. When D becomes bigger {D = 3) entanglement 
decreases slowly from C12 = 1 (Fig. 5d). 

We compare the concurrence versus magnetic field for different temperatures, when 
D = (Fig. 6) and when D — 2 (Fig. 7). In both cases at T = the entanglement 
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c a) D=0.1 c b) D=1.118 c c) D=1.19 




Figure 5: Concurrence in XX Z model versus temperature for S = 2 and a) 
D = 0.1, h)D = 1.118, c)D = 1.19, d)D = 3 

vanishes abruptly as B crosses critical value Be = V B^ + J-^ + Jz. This special point 
T = 0, B = Be at which entanglement becomes nonanalytic function of B, is the point 
of quantum phase transition. Comparison of figures 6 and 7 shows that the critical 
value of B at which entanglement disappears suddenly is growing with increasing 
couphng D: in Fig. 6, Be = 2 and in Fig.7, Be = 3.3. It shows again that increasing 
DM coupling improves entanglement. 



C 




Figure 6: Concurrence versus magnetic field B ioi D = and T = 0.1, 0.5, 1. 



8 XYZ Heisenberg Model 

In the present paper we are not going to analyze all possibilities for the XYZ model. 
Instead we restrict ourselves with a particular range of parameters to study the influ- 
ence of DM coupling in details. 
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Figure 7: Concurrence versus magnetic field B iov D = 2 and T = 0.1, 0.5, 1. 



8.1 Antiferromagnetic Case 
8.1.1 Pure Antiferromagnetic Model 

We start from the pure XYZ model, where for determinacy we chose Jz > Jy > Jx > 
implying J+ > |J_| > 0, J- = — jJ-l < 0. Eigenstates of the Hamiltonian ((2]) are 
Ei^2 ~ i l-^-l Mid E-s,4 = — =F |«^+|- For zero temperature the ground state is 
maximally entangled Bell state |01) — jlO) with the energy Eg = — \J+\- When 

temperature increases, the state mixes with higher states decreasing entanglement. To 
find concurrence we have the biggest eigenvalue X4 = ^ exp ^ "^^ ^2kT^'*' ^^'^ 

„ sinh — cosh ■^e""'^^''*^ , , , 

C12 = max{ ^ ,0}. (52) 

cosh -jji + cosh xt^^''^''*'^ 

Then entanglement occurs when sinh ^ > cosh ■^e~''^'"°^. It shows that entangle- 
ment depends essentially on the anisotropy, and grows with J-|- and decreases with J_ 
|28l. 



8.1.2 XYZ Model with DM Coupling 

Inclusion of the DM coupling, remains the energy levels Ei and E2 the same as above, 
while E3,4 = --1^ T In this case the ground state continues to be entan- 

gled state but with the energy E3. With growing temperature, mixing of this state 
with the higher states decreases the entanglement. If we consider the difference be- 



tween two lower states E4, — Eg = y + D'^, then by increasing the coupling D, it 
can be made arbitrary large, so that the entanglement will increase. For D >> \J+\ 
the state would be maximally entangled. At the concurrence T = Q 



C12 = { 



1, y^L»2 + J2 > J_ - J, 

1, Jd-^ + jI <J--J, 



(53) 



is nonanalytic function in D, and it signals about the quantum phase transition at 
D — Dc where 1 / + J\ — J- — Jz- When the temperature increases, entanglement 
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occurs for 

sinh > e-'-'"^ cosh ^, (54) 

and the concurrence 

^ sinh - e'-''^*''^ cosh , ^ 

Ci2 = (55) 

cosh + e--'*/*^'^ cosh -^jr 

increases with growing anisotropy J+ and the coupUng D. 

8.2 Ferromagnetic Case (J^ < < < 0) 

8.2.1 Pure XFZ Model 

Let Jz < Jy < Jx <Q then J+ = -| J+|, J- = | J- 1 > and = -\Jz\- For pure XFZ 
model, eigenstates of the Hamiltonian are E\^2 = — T l-^-l and £^3,4 = ± | J+|. 
For zero temperature the ground state is maximally entangled Bell state ]00) — 
with the energy Ei = — — | J_ | . With increasing temperature this state mixes with 
other states and entanglement decreases so that the concurrence 

sinh cosh ^e-l^'l/'=^ , ^ 

C12 = . 56 

cosh 1^ + cosh i^e-l-'^IAT- 

When temperature reaches the critical value T = Tc, given by a solution of the follow- 
ing transcendental equation 

sinh ^= cosh l^e-l^^l/'=^, (57) 

the concurrence vanishes and state becomes unentangled. 

8.2.2 XYZ Model with DM Coupling 

With inclusion of the DM coupling, the first couple of energy levels is the same i5i,2 = 
=F I J- 1 while the second couple becomes E3,i = T \J J+ + -D^- For D < Dc 

where Dc satisfies the equation D'i + = \Jz\ + \J-\, the ground state of the system 
is the maximally entangled Bell state |00) — If wo increase D, the difference 

between energy levels Ei and E3 decreases, so that a.t D = Dc the ground state 
becomes degenerate and entanglement vanishes. When D > Dc the ground state E3 
becomes entangled again. 

Due to the mixture of states by increasing temperature the entanglement decreases, 
so that, in the under critical region D < Dc the concurrence is 

sinh cosh ^^S+^e-l^^l/'^^ , , , 

C7i2 = max{ ^2 ,0}, (58) 

cosh %i + cosh ^^i—yp — e-l''»l/'=T 

k 1 k1 

while in the over critical region D > Dc it is 

sinh cosh ^ ^ ^ 

C12 = max{ ^, 0}. (59) 

cosh V ++ + e\JMkT cosh 1^ 
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For D — Dc, due to Ai = A3, the entanglement vanishes for any temperature. 
The entanglement dependence on T and D is shown in Figs. 8 and 9. For T — the 
figures show nonanalyticity at D = Dc which signals a quantum phase transition. 
The entanglement behavior in the under and the over critical regions is qualitatively 
different. For the under critical case with fixed temperature the entanglement decreases 
with growing D, and the level of entanglement quickly decreases with temperature. 
From another side, for fixed temperature in the over critical region the entanglement 
increases, and the level of entanglement decreases with temperature quite slowly. In 
addition if at T = we have only one critical point D = _Dc in which entanglement 
is zero, for T > entanglement vanishes at some interval which includes Dc and this 
interval extends with growing temperature. This is a result of ground state mixture 
with higher states. However by increasing D we can always lower the level of our 
ground state to decrease this mixture and increase entanglement. 



c 




Figure 8: Concurrence in ferromagnetic XYZ model versus coupling D at tem- 
perature T = 0.1, 0.5, 1 



9 Conclusion 

In conlusion we like to stress several implications for future studies. As was shown 
by Alcaraz et al. (35) the XXZ quantum Heisenberg chain with the DM interaction 
is equivalent to the pure XXZ Hamiltonian with modified boundary conditions and 
the anisotropy parameter, so that with these boundary conditions the model is still 
solvable by the Bethe Ansatz. Taking into account our results it shows possibility to 
control entanglement in XXZ model by varying boundary conditions. 

Recently it was found that the DM interaction can excite the entanglement and 
teleportation fidelity by using two independent Heisenberg XXX chains [36]. More- 
over, studying the effect of a phase shift on amount transferable two-spin entanglement 
[37] . it was shown that maximum attainable entanglement is enhanced by the DM in- 
teraction. Very recently geometric computations for a spin chain model with the DM 
interaction has been discussed in [32]. Finally it was found that the DM interac- 
tion is present in number of quasi-one dimensional magnets and is dominating for the 
compound RbCoCl3.2H20. These indicate that DM interaction could be significant 
in designing the spin-based realistic quantum computers [351. The above mentioned 
results suggest to study the most general XYZ Heisenberg model with DM interac- 
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Figure 9: 3D plot of concurrence in ferromagnetic XYZ model versus coupling 
D and temperature T 



tion as a quantum channel for quantum teleportation. These questions are now under 
investigation. 
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